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ABSTRACT 

A mathematical model linking a solute transport simulation 
model and a chance constrained, mixed integer programming 
optimization model for designing an optimal groundwater quality 
monitoring network is developed. 

The optimization model explicitly considers uncertainties of 
transport simulation and specified reliabilities of predicting the 
actual concentrations in time and space. The model is evaluated 
for different degrees of uncertainties in the transport modeling 
process by introducing cumulative distribution function <CDF > of 
the actual concentrations as constraints in the optimization 
model. The nonlinearities due to incorporation of CDF's are 
removed by using piecewise linearization scheme. The developed 
model is solved for various degrees of uncertainties and -for 
different reliability values. The solution results clearly show 
that the optimal well locations are dependent on the degree of 
uncertainty and the specified reliability values in the prediction 


af actual concentrations. 
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CHAPTER 1 

INTRODUCTION 

1.1 MOTIVATION FRAMEWORK AND OBJECTIVES 

There has been an increasing concern throughout the past two 
decades to protect our groundwater resources from contamination. 
The water supplied for domestic and agricultural uses may be 
polluted by any one of the following potential groundwater 
contamination sources, 

1) leaky effluent pipe lines, 

2) toxic chemical waste disposal sites, 

3) waste water infiltration, 

4) irrigation run-off infiltration, 

5) underground fuel storage tanks, and 
unidentified geologic formations. 

Identification of groundwater pollution source and monitoring 
of groundwater quality are essential steps in groundwater quality 
management. In this study a mathematical model for designing a 
groundwater quality monitoring network is deve loped that links a 
groundwater pollution transport simulation model and an 
optimization model, explicitly considering uncertainties of 
transport simulation. To illustrate the performance of this 
model, tritium is considered as the (radioactive) pollutant. 
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Nuclear reactors use heavy water as moderator and coolant 
which results in considerable production of tritium as an 
activation product. During their normal operation tritium is 
routinely released into the environment from the reactors through 
atmospheric and liquid discharge routes, generally within 
tolerable discharge limits. On release to the environment, 
tritium <T ^ = 12.3 years) being an isotope of hydrogen 
assimilates readily with water component of the atmospheric, 
aquatic and biological systems. In case of an accident large 
quantity of tritium may enter the groundwater system and can 
impair the use of water. 

In contrast with surface water pollution, subsurface 
pollution is difficult to detect, is even more difficult to 
control, and may persist for decades. Thus with the growing 
recognition of the importance of groundwater resources, efforts 
are increasing to prevent, reduce, and eliminate groundwater 
pollution. It becomes essential for the water manager to take 
this responsibility of detecting the harmful contaminants which 
moves very rapidly in the saturated zone. Concentration contours 
at different time periods for pollutant disposal from two sites 
is shown in Figure 1—4. Remedial actions to combat the threat of 
contamination requires detection of pollution sources and spatial 
extent of such concentration. It is therefore necessary to design 
a groundwater quality monitoring network so as to detect and 
contain the extent of contamination. 
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A chance constrained, mixed integer programming model is 
developed for designing an optimal groundwater quality monitoring 
network, incorporating uncertainties in the prediction of 
pollutant movement in the saturated zone. The nonlinearities due 
to the incorporation of cumulative distribution functions are 
accommodated through a piecewise linearization scheme. The design 
of this optimal monitoring network is based on the solutions of 
two mathematical models ; a simulation model for prediction of 
salute transport in the saturated zone, and an optimization 
model. The optimization model utilizes the solutions from the 
simulation model. In particular, the constraints of the 
optimization model are constructed by incorporat ing simulation 
results. The simulation model provides information about the 
pollutant transport with respect to time and space. The chance 
constrained optimization model specifies the optimal location of 
the monitoring wells subject to the maximum limit on the number 
of such wells. 

Parameter estimation uncertainties, heterogeneity of the 
porous media etc. need to be incorporated while simulating 
groundwater contamination due to known sources. Optimal design of 
a monitoring network is a difficult task because, it is difficult 
to correctly predict the movement of pollutant in groundwater. 
Also, using large number of groundwater quality monitoring wells 
is not always physically and/or economically feasible. Therefore 
there is a necessity of optimal decisions regarding the design of 
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a monitoring network. These optimal designs should also consider 
inherent uncertainties in the prediction of pollutant transport. 

The primary aim of this work is to design an optimal 
groundwater quality monitoring network incorporating parameter 
estimation errors and related inaccuracies in the modeling of the 
transport process. These errors and uncertainties in modeling may 
result from errors in observed groundwater table or hydraulic 
head. These errors and uncertainties in simulation and prediction 
of salute transport in groundwater can be indirectly accounted for 
by randomly varying the response matrix, representing influence of 
pollutant sources upon spatial concentration values evolving over 
time. The link between the simulation and optimization model is 
provided through the concentration response matrix. 

The role of chance constraints in the optimization model is 
to introduce a measure of reliability in the predicted 
concentrations at specified location. Increasing this reliability 
will result in more conservative prediction of the spatial and 
temporal distribution of the resulting pollutant concentrations 
due to known sources. The chance constraints are based on 
specified values of these reliabilities and the cumulative 
distribution function <CDF> of the actual value given an 
estimated value obtained by simulation. The CDF are incorporated 
in the optimization model. The CDF's of the spatial and temporal 
values of the concentrat ions can be obtained by random variations 
in the estimated inputs to the simulation model, or by random 
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variations of the response matrix that represents the response of 
the subsurface saturated zone for given input value of parameters. 
The general objective of the chance constrained optimization model 
is to maximize the probability of detecting contamination at 
locations, where the standard is exceeded, with due weightage 
given to the degree of exceedences. These exceedence values are 
computed in terms of the predicted values of concentrations 
expressed as a function of the specified rel iabi 1 i t ies. 

1.2 SPECIFIC OBJECTIVES OF STUDY 

Ci) Use of USSS (Konikaw and Bredehoeft, 1978) computer 

model for solute— transport simulation. 

<ii) Formulate a chance constrained optimization model 

linked to a simulation model, for groundwater 
pollution monitoring network design taking into 
account uncertainties in the estimation of 
hydraulic and transport parameters, and other 
modeling errors. 

(iii) Test the performance of the optimization model for 

different degrees of parameter uncertainties for 
specified rel iabi 1 i t ies. 

1.3 BACKGROUND AND SURVEY OF PREVIOUS WORK 

An excellent review of models for groundwater quantity and 
quality management is given in the classic paper by Borelick 
(1983). Much of the following discussion is based on this work 
with appropriate additions wherever necessary. 
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i.3.1 Groundwater- Modeling 

With the help of numerical groundwater flow and transport 
modeling, it has now become much easier to understand the 
groundwater systems. Simulation of a groundwater system refers to 
the formulation and operation of a model whose behaviour assumes 
the appearance of the actual aquifer behaviour. According to 
Mercer and Faust <1980), groundwater models, dealing with 
groundwater flew hydraulics and solute transport can be placed 
under four general classifications. The models and their 
application are listed in Figure 5. 

In order to take care of uncertainties in the estimation of 
model parameters the' groundwater models can be combined with 
statistical techniques. Data preparation for the groundwater 
model first involves determining the boundaries of the model to be 
modeled, thereafter aquifer parameters and initial data for each 
grid such as dispersivi ty , hydraulic conductivity, porosity are to 
be specified. Care should be taken while using a numerical model 
as these models are based an a set of simplifying assumptions 
which limit their use for specific problems only. 

i . 3. 2 Groundwater Management 

In the past, numerical simulation models have been used to 
evaluate groundwater resources. Numerical simulation models that 
solve groundwater flow or solute transport equation in conjunction 
with optimization techniques like linear or quadratic programming 
are powerful aquifer management tools. 
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Applications 


• WATER SUPPLY • SEA-WATER • GEOTHERMAL • LAND SUBSIDENCE 

INTRUSION 

• REGIONAL AQUIFER • LAND FILLS • THERMAL 

ANALYSIS STORAGE 

•NEAR-WELL • WASTE • HEAT PUMP 

PERFORMANCE INJECTION 

• GROUND-WATER/ • RADIOACTIVE • THERMAL POLLUTION 
SURFACE WATER WASTE STORAGE 

INTERACTIONS 

• DEWATERING • HOLDING PONDS 

OPERATIONS 

• GROUND WATER 
POLLUTION 

Figure 5". Types ot ground-water models and typical 
applications^ Mercer and Faust, .1980) 
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Simulation as both a method to explore hydrogeolog ic 
problems and a tool to predict impacts upon groundwater system 
will continue to be essential to hydrologists and to water 
managers. Simulation models are often utilized to explore 
groundwater management alternatives. In such cases a model is 
executed repeatedly under various design scenarios which attempt 
to achieve a particular abjective, such as preventing saltwater 
intrusion, dewatering an excavation area or isolating a plume of 
contaminated groundwater. In such an approach often physical and 
operational restrictions are not considered in groundwater 
management goals. So a joint simulation and management model 
approach is needed. 

Borelick <1982> classified groundwater hydraulic management 
models into two basic types 

i) Embedding method 

ii) Response Matrix approach. 

The embedding method for the hydraulic management of aquifer uses 
linear programming formulation that incorporates numerical 
approximations of the groundwater equations as constraints. In the 
response matrix approach an external groundwater simulation model 
is used to develop unit responses. The assemblage of unit 
responses form the response matrix is included in the management 
model. In the embedding approach the flow equation are included 
in the linear programme as constraints and a complete simulation 
is solved as a part of the optimization model. This helps in 
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getting a great deal of aquifer information. But it is not always 
desired to involve all the hydraulic heads over time and space. 
Hence many of the decision variables and constraints are not 
necessary in the linear programming or other optimization models. 

In response matrix approach solutions to the flow equations 
serve as constraints. It is very efficient as compared to 
embedding method as constraints are included only for specified 
locations and times. So the response matrix approach can be used 
for large transient systems in an efficient manner. 

1.3.3 Groundwater Quality Management 

The primary aim of the groundwater quality management model 
is to minimize the harmful effects of waste disposals by 
maintaining the water quality standards. The joint use of 
numerical simulation and optimization models has been applied for 
groundwater pollutant source management. The need here is to 
utilize the aquifer for both waste disposal and for water supply. 
In the fallowing section steady state groundwater quality 
management model and transient groundwater quality management 
model are discussed. 

1.3. 3. A Steady state groundwater quality management model 

In steady state groundwater quality management model, the 
sources of pollution do not change with time, so the problem is 
only space dependent. Here the storage coefficient is set equal 
to zero. The time derivatives in the governing equations are also 
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set equal to zero. In groundwater aquifers steady state 
conditions normally do not exist. 

1.3. 3. B Tr ansient groundwater quality management model 

In case of transient groundwater quality management the 
pollution sources occur over both space and time. So the 
identification of the pollutant source and the time over which the 
pollution occurs are equally important. The concentration 
measurements are collected over time at various locations. 
Although the transient case models are more computationally 
complex they tend to represent more accurately the aquifer being 
studied. 

1.3.4 Response Matrix Approach in Groundwater Quality Management 
Problems. 

By definition, the concentration response matrix ER3 
describes the influence of a unit disposal flux on the 
observation well concentration as a result of the unit pollutant 
injections at the potential source sites. A groundwater solute 
transport model is used to simulate these response concentrat ions 
for the CR3 matrix. These simulation results produce a set of 
breakthrough curves for each potential disposal site. Each set 
contains curves for each of the observation wells mampping the 
concentration resulting from disposal at a particular site. 

These curves resemble a unit hydrograph in shape. In 
actuality, the principle upon which they are based is same. As 
with unit hydrograph theory, a unit disposal, flux is applied to 
the groundwater system to simulate the resulting response curve. 
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Each set of curves is then referred to as a suite of breakthorough 
curves which are used to construct the concentration response 
matrix CR3. Some of these breakthrough curves are shown in 
F igureS. 

1.4 GROUNDWATER QUALITY MONITORING NETWORK DESIGN 

In order to prove that the groundwater is contaminated or to 
formulate effective remedial measures, it is essential to test the 
water samples from the observation wells that are located in the 
affected area. However economic and other constraints limit the 
number of wells that can be installed and it becomes essential to 
design an optimal groundwater quality monitoring network. 

With the help of groundwater salute transport simulation 
model, a likely transport scenario can be obtained. The 
uncertainties in the estimated parameter values, and the 
observation errors in hydraulic heads are not incorporated in the 
simulation model. Moreover for economic feasibility it is not 
passible to construct a large number of observation wells. It is 
necessary to apply an optimum criteria to the selection of a 
monitoring network design. An optimization model can be used to 
solve these problems. Such an optimization model must consider 
all physical and managerial constraints and have an explicit 
abjective function. 

An approach to the design of a groundwater quality monitoring 
network was suggested by Meyer and Brill <1988). This method can 
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observation: well locations 
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be used to select a network that maximize probability of detection 
in the -face of uncertainty. While evaluating the performance of 
the model as proposed by Meyer and Brill (1988) a number of 
shortcomings were detected, and a new model was suggested by 
Datta arid Purwar 11991). A reverse approach of utilizing the 
probability of detection for a given monitoring network was 
suggested by Massmann and Freeze (1987), as part of a larger 
problem for design of a landfill operation. 

The model proposed in this study is capable of designing an 
optimal groundwater quality monitoring network, under condition of 
uncertainty and incorporates measures of reliability in the 
prediction of actual concentration at a given location and time. 
This model is formulated as a mixed integer programming with 
chance constraints. 

Chance constrained model, a particular form of stochastic 
optirai zation model has been used in various areas of water 
resources management. Houck (1979) developed a chance constrained 
model fear reservoir management. Tung (1986) presented a chance 
constrained model for groundwater quantity management. This model 
explicitly considers the random nature of transmissivity and 
storage coefficient, that enable the determination of optimal 
pumping pattern in a well field subject to a specified system 
performance reliability. This work also utilizes first order 
uncertainty analysis for hydraulic conductivity and storage 
coefficient. However, this model is not applicable to groundwater 


quality management. 
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CHAPTER 2 

tCTHODOLOGY 

2.1 DESCRIPTION OF THE METHODOLOGY 

The proposed chance constrained optimization approach to 
design an optimal groundwater quality monitoring network under 
conditions of uncertainty has two major components. The first 
component involves the use of a groundwater solute transport 
computer simulation model to simulate pollutant movement from 
known sources to potential groundwater quality monitoring wells. 
In the second component simulated concentration data incorporating 
specified measures of uncertainty in the modeling process are 
utilised to formulate the chance constrained optimization model 
for solving a groundwater quality monitoring network design 
problem. 

The proposed model includes a nonlinear abjective function 
and some nonlinear constraints due to cumulative distribution 
function that are introduced explicitly in the model. However, an 
equivalent formulation of this nonlinear optimization model is 
presented and used in this study, using piecewise linearization 
technique and mixed integer programming. The groundwater solute 
transport simulation model, the optimization model and the 
procedure for constructing the response matrix is discussed in 
detail in the following sections. 
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2.1.1 USGS Two-Dimensional Flow and Solute Transport Model: 

The solute transport in groundwater was simulated using a 
finite difference based 2-D, numerical model by Konikow and 
Bredehoeft (1978). The USGS computer model for solute transport 
combines the groundwater flow equation with the solute transport 
equation in order to compute the transient changes in 
reactive/nonreactive solute concentration. The model may be 
applied to both steady state and/or transient flow problems for 
one dimensional or two dimensional flows. Its purpose is to 
identify the change in dissolved chemical concentration in the 
aquifer at desired locations over time. Due to the process of 
convection, hydrodynamic dispersion, chemical reactions, 

adsorption, fluid sources or sinks, and dilution or mixing of the 
fluid sources, there is change in concentration . 

The hydraulic head distribution in the aquifer is simulated 
by the groundwater flow equation. Pinder and Bredehoeft (1968) 
derived the following equation to describe the transient two 
dimensional areal flow through a non-homogeneous , anisotropic 
confined aquifer by a homogeneous incompressible fluid (using 
Einstein notation) 

( T ii £t] = 5 & + w = ‘* 2 u> 

x J J 

2 

where* T. . = transmissivity tensor* L /T ; 
x Jf 

h = hydraulic head, L; 

S = Storage coefficient, (dimensionless); 

t = time, T; 
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W = W (x,y,t), volume flux per unit area (positive far 
outflow and negative for inflow), L/T; and 
x , x = Cartesian coordinates, L. 

1 J 

If the fluxes, W<x,y,t> deal with direct withdrawl or recharge or 
a steady leakage into or out of the aquifer, then the following 
equation is used: 


k 

W <x,y,t> = Q(x,y,t) — < H — h ) (2) 

m s 

where, Q = rate of withdrawl (+ sign) or recharge <- sign), L/T; 

k^ = vertical hydraulic conductivity of the confining 
layer, streambed, or lakebed L/T, 
m = thickness of the confining layer, streambed or 
lakebed, L; and 

= hydraulic head in the source bed, stream, or lake, L. 
The derivation from Darcy's law for the average seepage velocity 
is as fallows. 


V. 

1 


k. . . 

1 j <9h 

e &.< . 

3 


(3) 


where, = Seepage velocity in x direction, L/T; 


k_ = hydraulic conductivity tensor, L/T; and 
e = effective porosity of the aquifer (dimensionless) 

The chemical concentration in the aquifer system is described 
by solute transport equation. The following is the equation for 
transient two-dimensional areal transport and dispersion of a 
given single reactive dissolved chemical species in the flowing 
groundwater CSorelick, et al., 19833 

sc <? r _ sc > 


±1 * I D 1 

£t . ! i j Sx . J 

1 v J 


* C'W 

(C V. ) - 

Cx . 1 eb 

x 


>X, i, j = 1,2 (4) 
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where, 

C = Concentration of the dissolved chemical species; mg / 1 
P ci/t; 

D- • = dispersion tensor, m /d 

* J 

= average pore water velocity in the direction x^ , m/d; 

b = saturated aquifer thickness; m; 

e = effective aquifer porosity; dimensionless; 

C' = solute concentration in fluid source or sink, mg/L or 
p ci/1; 

W = volume flux (source) per unit area, m/d; 

X = first order kinetic decay rate, 1/d; 

x., x. = cartesian coordinates, m; 

13 

t = time, d. 

The first term on the right hand side of equation (4) describes 

the change in concentration due to hydrodynamic dispersion. The 

effects of convective transport are shown in the second term on 

the right hand side and fluid sources and/or sinks are represented 

by the third term the fourth term shows the effect of radioactive 

decay on the change in concentration . According to Bear (1972), 

hydrodynamic dispersion is the macroscopic outcome of actual 

movements of individual tracer particles through the pares. This 

includes the process of mechanical dispersion which is dependent 

on fluid flow through the system and the process of molecular and 

ionic diffusion which are independent of flow. In the development 

of the USGS model, the definable contribution of molecular and 

ionic diffusion to hydrodynamic dispersion is assumed negligible. 

The dispersion coefficient (D. .) is then defined using 

1 3 
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Sche idegger' s C19613 equation relating the coefficient to the 
velocity of ground water flow and to the nature of the aquifers 


a. . 

1 jm,n 


V V 
m n 


where : 


i jmn 


= dispersivity of the aquifer; L; 


V s components of velocity in m and n directions 
respect ively , L/T and 

jV{ = the magnitude of the velocity, L/T. 

In addition, Scheidegger (1961) defines the dispersivity tensor in 
terms of two constants for an isotropic aquifer. The longitudinal 
and transverse dispersivit ies of the aquifer (a^ and a^. 
respectively) are related to the longitudinal and transverse 
dispersion coefficients and D^> as follows z 

D l = a L [V{ (6> 

and D T = a T {V{ (71 

There are numerous assumptions associated with the use of the flow 
and solute transport equations. These assumptions must be taken 
into account when applying the USGS solute transport model to an 
actual field problem. The assumptions summarized by Konikow and 


Bredehoeft Cl 978 3 are as follows: 


1. Darcy's law is valid and the only significant mechanism for 
driving the flow through the system are the hydraulic— head 
gradients. 

2. Porosity and hydraulic conductivity of the aquifer are 
constant with time and porosity is uniform with space. 

3. The velocity distribution is not affected by the gradients of 
fluid density, viscosity, and temperature. 
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4. Vertical variations in head and concentration are negligible. 
5- The aquifer is homogeneous and isotropic with respect to the 
coefficients of longitudinal and transverse dispersivity. 

6. Ionic and molecular diffusion contribution to the total 
dispersive flux are negligible. 

The computer model developed by Konikow and Bredehoeft C19783 is 
based on a rectangular, block— centered, finite difference grid 
with nodes at the center. This allows for the input of saturated 
thickness, transmissivity, boundary conditions initial heads and 
concentrations, observations points, injection or withdrawl wells, 
and spatially varying diffuse recharge or discharge conditions. 
Because analytical solution to the flow and salute transport 
equation's cannot be salved directly due to the variable 
properties and complex boundary conditions of an aquifer, the 
model uses an alternating direction implicit (ADI) method to solve 
the finite difference equations describing the groundwater flow 
equation and the method of characteristics to solve the solute 
transport equation. 

2.1.2 The Concentration Response Matrix 

To avoid the repeatitive use of simulation model and to 
provide the necessary input to the optimization model the response 
matrix (Gorelick, 1982) approach has been used. This particular 
approach is discussed here. The concentration response matrix 
describes the influence of unit pollutant sources upon 
concentrations at specified locations over time. 

Each simulation far a unit injection rate at the source 
results in a suite of breakthrough curves at the groundwater 



22 


quality monitoring sites. To create other realistic or probable 
concentration response curves random error terms are added to each 
element of the response matrix to account for errors in estimation 
of groundwater field parameters, and other modeling uncertainties. 

Each suite of breakthrough curves serves as one column of the 
concent rat ion response matrix. A complete assemblage of these 
constitutes the concentration response submatrix for one 
management period. This single management period response 
submatrix is called R^. Response submatrices representing 
additional management period's may be developed by staircasing the 
original response submatrix. The original submatrix R t is 
repeated and shifted down to account far successive periods. The 
concentration response matrix CR3 therefore consists of several 
submatrices R^, R^, . - R p where p represents last management 
period. The general structure of CR3 matrix is displayed in Fig. 
7. The time varying pollutant concentrations at the monitoring 
sites for different injections of pollutants with known 
concentrations at different management periods are given by 

CR3 Cf 3 = CC3 <8> 

where 

CR3 the response matrix 

Cf 3 the known actual disposal at sources 

CC3 concentrations at monitoring wells 1 .. n 
Therefore once the CR3 matrix has been computed using a 
mathematical groundwater transport simulation model, the resulting 
concentrations CC3 at different times and locations can be 
obtained using equation <B) for a given Cf3. This eliminates the 




Figure 7. Structure of a concentration response 
matrix for a ground water quality 
monitoring wells. 




necessity of repeatedly using the simulation model to compute CC3 
for any new pollutant input condition represented by CfJ. 

2.1.3 Computation of the Response Matrix 

Each simulation of a unit injection rate resulted in a suite 
of breakthrough curves at the monitoring wells. Each suite of 
breakthrough curves serves as one column of the concentration 
response matrix. These breakthrough curves (Gorelick, 1982) for a 
few selected potential observation well locations are shown in 
Figure £>. These breakthrough curves were constructed for each of 
the ten potential groundwater quality monitoring well sites for 
each of the two source locations. 

In this study, two pollution sources are considered, 
therefore, the complete collection of two columns (one for each 
potential disposal site and S^) constitute a submatrix 
containing the concentration response data for 8 years simulations 
period at each of ten potential observation well location, 1 
through 10. the single management period for constructing the 
response matrix and for the optimization model is assumed to be of 
1 year duration. Each submatrix is shifted down by one management 
period. 


2.1.4 Incorporation of Uncertainties in Solute; Transport - 

Simulation i 

The response matrix can be utilized for simulating the 
pollutant migration without using an analytical or numerical 
simulation model to predict resulting concentrtion for different 
sources or input conditions. However, these simulated transport 
conditions are based only on deterministic conditions, ignoring 
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the uncertainties in parameter estimations and modeling errors. 

These uncertainties are incorporated here by random perturbat ion 

of the CR3 matrix, using statistically generated error terms. 

Thus the originaally computed matrix is assumed to be uncertain 

and then it is varied by adding randomly generated variables to 

each element of the matrix. This is done by adding random number 

generated with specified means and variances. A random error 

term, £ is added to each elements of the response matrix CR3 to 

obtain new values for the response matrix CR'3 

CR' 3 — ER3 -i- 3 <9> 

and CR # 3 Cf3 = CC3 <10> 

The random error terms given by each element of C£3 is defined as 

. 3 = e a + v (11 ) 

* •• j 

where the generted random variable Cel is sampled from a Normal 
Distribution with a mean p equal to zero and standard deviation O', 
equal to a x R. where a is a fraction. 

2. J 

for ot = 0.25, and = 0.0 equation <111 becomes 

R' . = R. . + e. . * (0.25 R. . > (12) 

u u u n 

A FORTRAN programme for statisticl perturbtion of the CR3 matrix 
is given in Appendix A. 

2.1.5 Computation of the Cumulative Distribution Function of 
Actual Concentration 

The cumulative distribution functions (CDF's) of actual 
concentrat ions are computed using statistically generated values 
of CR'3 and utilizing equation 10. For each given CR'l a value of 
probable concentration at a specified location is generated. For 
a large set of generated CR'3 matrix, a large number of 
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concentrations are generated statistically. These CO values far 
a particular location at a specific time period form a sample set 
of concentration at a location and time, for the same input 
condition. 

The plotting position formulae (or other appropriate methods) 
can be used to obtain the CDF for actual concentration at a 
location at a given time. This method is adopted in this study. A 
typical cumulative distribution function for concentrations at a 
given location for different degree of uncertainty (a) is shown in 
Figure 8. 

2.1.6 Chance Constrained Optimization Model for Groundwater 
duality Monitoring Network Design 

A chance constrained optimization model is presented in this 
section for an optimal design of groundwater qualiaty monitoring 
network. The role of the chance constraints in the model is to 
introduce a measure of reliability in the estimated values of 
predicted concentrations at specified locations. The estimated 
values of predicted concentration depend upon the degree of 
uncertainty incorporated. 

The optimization model incorporates the cumulative 
distribution function (CDF) of the actual concentration for each 
potential groundwater quality monitoring well location for 
specified management periods. The CDF's of the actual 
concentrations can be obtaianed by randomly varying the response 
matrix, that represents the unit response of the subsurface 
saturated zone for given input conditions. 




Concentrations 

Figure 8 Typical cumulative distribution function 


for concentrations at a given location 
for different values of a 
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The general abjective of the chance constrained optimization 
model is to maximize the probability of detecting contaminations 
at locations, where a given polution standard is exceeded. This 
general objective is mathematically stated in this model as 
minimization of locationwise sum of the positive difference 
between the undetected spatial and temporal pollutant 
concentrations and the contaminant concentration standard, subject 
to imposed chance constraints. 

This nonlinear objective function is converted to linear 
objective function with suitable modifications in the constraints. 
The nonlinear constraints due to nonlinear CQF's introduced in the 
model are transformed to linear constraints by using piecewise 
linearization technique. The resulting mixed integer programming 
model is solved using LINDQ (Linear Interactive Discrete 
Optimizer ) . 

A chance constrained mixed integer programming model for the 
monitoring network design can be stated as: 

Minimize: Z = E E <x k - S .) (1 - y > (13) 

k«K'j«J' 3 3 3 


subject to : 

E y i 
j J 3 

< P 

(14) 

P cc k < 

x k 3 = ft k , for V j e J 

3 3 

(15) 



fl k £ r k : 

> 

* 

**■* 

> 

(16) 

J J 

X k > O 


(17) 


3 


(O, 1) 


(18) 
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T he 


fol lowing 

k 

;< . 

J 



S. 

J 


<E 





J,K 


natation are used. 

: probable concentration at potential monitoring site 

j at management time period k, expressed as a 

function of reliability ft:. 

J 

: actual concentration at site j at management period 
acceptable standard concentration at potential 
monitoring site j. 

: belonging to the set 

: a decision variable indicating if a monitoring well 
is to be installed at site j, a value of 1 
indicating installation. 

maximum number of wells permitted for the entire 
monitoring network 
for all values of 

inverse CDF for the predicted concentrations at 
site j at management time period k. 

probability that the actual concentration at site j 
at management time period k will not exceed a 
particular predicted value. 

a specified value of reliability denoting the 
required probability that a predicted value of 
concentration at site j at time period k will not 
be exceeded. 

the set of possible value of j & k. 

subsets of J and K containing those values of j %/. k 


k. 


J' ,K 


Jc 

for which <x. - S.) are positive. 
J J 
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In the above formulation the objective function minimizes the 
sum of the positive differences between predicted concentrations 
at potential monitoring sites and the imposed concentration limit 
or standard at all potential monitoring sites where a monitoring 
well is not to be installed according to the solution of the 
model . 

This objective function is nonlinear due to the product of 

k 

two decision variables :<.* and y. occurring in this function. Also 

w J 

it is not efficient to determine externally, all the sites where 
the concentration standard is exceeded at one or more management 
period . 

Constraint (14) limits the total number of monitoring wells 

to less than or equal to a permissible maximum value of P. 

Constraint set <15> defines the probability or the reliability 

that the actual concentration at site j and management period k 

will not exceed a particular value. Constraint set (16) relates 

these probabilities to the specified reliability value y.'. 

Constraint set (18) denotes that y. are (0,1) decision 

J 

variables, a solution value of y . equal to zero will mean that a 

«! 

monitoring well is not to be installed at site j , while a value of 
y^ equal to 1 will specify that a well is to be installed at site 

j - 

The chance constraint set (15) can be replaced by an 
equivalent constraint 

x k = F" 1 (/£> V j € J <!*?> 

j Ok j 

k 

Here x . is expressed as a. function of the probability that the 
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actual value of the concentration at site j and management period 
k will not exceed a certain predicted value* This value can be 
obtained from the CDF's of the actual concentrations, that are 
computed by considering the uncertainties in the estimation of 
ac tua 1 concen t r at ions . 

The above mentioned model is difficult to solve due to the 
nonlinearity of the abjective function and also due to the need to 
identify the subsets J' and * Also the solution is further 
complicated by the decision variable y., which can take only 
integer (O, 1) values. Therefore, this model is modified as 
follows to convert it into a linear, mixed integer programming 
model. The modified model can be stated as 

Minimize: Z = £ E D k (20 ) 

k <E K j € J 3 

subject to : 


D k > x k - S. - My. V j, V k 
j a j j 

E y * P 
j <= J J 

v j * 1 ' 


(21 > 

(22 > 

(23) 


P k r k 

i J 


(24 > 


k k 

D*, x* > 0 

a a 


y. = (O, 1) 


(25) 

(26) 


The additional notation used is, M = a positive large real number. 
Due to the introduction of constraint set (20) the original 
objective function (13) can be converted to an equivalent linear 
objective function that also ensures that only the positive 
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if 

deviations i.e* positive value of S. values are included in 

the objective function, at potential locations where a monitoring 
well is not to be installed* 

The product of two decision variables occurring in the 

objective function (13> does not appear in the new objective 

function (20) * the large value of M ensures that whenever y is 

j 

nonzero, D * can take values equal to zero* Also* when :<*** — S is 

j j 

negative, even if y . is equal to l, D k can take a value equal to 

•J J 

zero* The formulation is still nonlinear because the CDF's of 

k 

actual concentrations are nonlinear functions of ft To overcome 

3 

this problem, piecewise linearization scheme is used to linearize 
the CDF's. Figure 9 shows a typical CDF linearized by three 
straight lines. The piecewise linearized CDF's are introduced in 
the optimization model using suitable constraints. The resulting 
model can be salved as a mixed integer programming model. The 
piecewise linearization scheme is explained in details in the 
fallowing section. 

The piecewise linearized CDF's are represented in the 
optimization model with the following new constraints, assuming 
that three linearized pieces are used to linerize a particular 
CDF, in the range of interest. 

ft. < m ik x k + cf V j, V K <27> 

a * * 

ft k < m^ W K k + C| k V j, V K <28> 

-k . jk k ^ _ jk w . w 

« 2 X j + C 2 V V K 


( 29 ) 




and location 

Figure 9 Piece wise linearization scheme for 
the decision space def ined by the 
cumulative distribution function of 
the pollutant concentrations. 
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where 

j Jj- J |(- j 

m l ’ *"2 » ia 2 ' ~ CQe ^ * icients of the linear equation representing a 

linearized portions of a CDF for a given j and k. 

j Jc* j |c; J k 

C i C 3 = constant terms in the linearizing equations for 

given values of j and k. 

t* 

It is not required ta consider the values of tfi.l less than 

J 

60 or 7054 so that the decision space can be safely assumed to 
remain convex, guaranting an optimal solution. The objective 
function <20 ) together with the constraints (21) - (26) is the new 
optimization model for solving the monitoring network design 
problem. This model explicitly incorporates the CDF's of actual 
concentration at different locations and management periods. 
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CHAPTER 3 

SPECIFIC APPLICATION AND DISCUSSION OF RESULTS 

3.1 STUDY AREA ^ 

The study area comprises of a portion of an aquifer. This 
study area is of irregular geometry with impermeable boundaries on 
the left and right, and nan uniform steady flow from top to bottom 
of the <Figure 10). The finite difference grid applied to the 
system is shown in Figure 10. The two sources of pollution S x and 
S 7 as well as the potential observation well locations are shown 
in Figure 11. The coordinates of observation well locations are 
given in Table 1. The aquifer is assumed to be homogeneous and 
isotropic. The aquifer parameter values are given in Table 2. 

A radioactive pollutant (Tritium, H , = 12.3 years) is 

assumed to enter the groundwater system from sources and S^. 

The actual disposal of the pollutant takes place during the first 
4 management period only. The pollutant fluxes from sources 
and far four management periods are given in Table 3. In this 
study the maximum permissible concentration for tritium is taken 
as 3 idci/l for all locations. 

3.2 DISCUSSION OF RESULTS 

The proposed optimization model was solved for a number of 
scenarios, with different degrees of uncertainties, and with 
different values of P and the reliabilities. Table 4 shows the 
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Table 1 i Potential Groundwater Quality Monitoring 
Well Locations. 


Well No. Grid Coordinates <I,J) 


1 

<12, 

7) 

2 

<12, 

8) 

3 

<12, 

9) 

4 

<13, 

7) 

5 

<13, 

9) 

6 

<14, 

5) 

7 

<14, 

6 > 

e 

<14, 

7 > 

9 

<14, 

9) 

10 

<14, 

io> 


Table 2s Aquifer Parameter values. 


Initial Concentration 
Aquifer thickness 
Effective Porosity 
Longitudinal dispersivity 
Transverse dispersivity 
Hydraulic Conductivity 


0.0 jJ ci/ I 
30.5 m 
0.2 


7.625 m 


3.05 m 


13 m/day 



3 *? 


Table 3: Actual Disposal Fluxes. 



Site No. 

Actual Disposal Flux (fjci/s) 

Year 1 

SI 

0.01 


S2 

0.01 

Year 2 

SI 

0.01 


S2 

0.0 

Year 3 

SI 

0.005 


S2 

0.010 

Year 4 

SI 

0.0 


S2 

0.005 
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Table 4: Solution Results of the Optimization Model. 


a 


0. 1 


O. I 


0.1 


0.1 


P Objective Optimal solution in 

function value terms of well 

location 


0.75 3 1 

2 

7 

•w-f 

4 

5 

6 

7 

8 
9 

lO 


202.9 

160.09 
118.90 

82.49 

47.52 

19.519 

0.5247 

0 

0 

O 


0.8 3 1 

2 

3 

4 

5 

6 

7 

8 
9 

lO 


207.63 

164.18 

121.41 

84.2 

48.52 

19.876 

0.6075 

0 

0 

O 


0.9 3 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


219.75 

173.56 

127.42 

88.87 

57.03 

20.729 

0.7729 

O 

O 

O 


0.95 3 1 

2 

3 

4 

5 

6 

7 

8 
9 


227.19 

178.96 

131.201 

91.81 

52.75 

21.57 

1.0 

O 

0 


lO 0 


t 

1,4 

1.4.9 

1.4. 8. 9 

1.3. 4. 8. 9 

1.3. 4.5. 8. 9 

1 . 2 . 3 . 4 . 5 . 8. 9 

1 . 2 . 3 . 4 . 5 . 7 . 8. 9 
1 , 2 , 3 , 4 , 5 , 7 , 8, 9 
1 , 2 , 3 , 4 , 5 , 7 , 8, 9 


1,4 

1 . 4.9 

1 . 4 . 8. 9 

1 . 3 . 4 . 8. 9 

1 . 3 . 4 . 5 . 8. 9 

1 . 2 . 3 . 4 . 5 . 8. 9 

1 . 2 . 3 . 4 . 5 . 7 . 8. 9 
1 , 2 , 3 , 4 , 5 , 7 , 8, 9 

1 . 2 . 3 . 4 . 5 . 7 . 8. 9 

t 

1,4 

1 . 4.9 

1 . 4 . 8. 9 

1 . 3 . 4 . 8. 9 

1 . 3 . 4 . 5 . 8. 9 

1.2. 3. 4. 5. 8. 9 

1 . 2 . 3 . 4 . 5 . 7 . 8. 9 
1 , 2 , 3 , 4 , 5 , 7 , 8 , 9,10 
1 , 2 , 3 , 4 , 5 , 7 , 8 , 9,10 

1 

1 . 4.9 

1 . 4 . 8. 9 

1 . 3 . 4 . 8. 9 

1 . 3 . 4 . 5 . 8. 9 

1 . 2 . 3 . 4 . 5 . 8. 9 

1 . 2 . 3 . 4 . 5 . 7 . 8.9 
1 , 2 , 3 , 4 , 5 , 7 , 8 , 9,10 
1 , 2 , 3 , 4 , 5 , 7 , 8 , 9,10 
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Table 4 (Continued). 


a 


Objective Optimal solution in 

function value terms of well 

location 


0.25 0.75 3 1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


239.39 

4 

187.07 

1,4 

142.58 

1,3,4 

98.34 

1,3, 4, 8 

60.07 

1,3, 4,8, 9 

24.07 

1,3, 4, 5, 8, 9 

1.87 

1,2,3, 4, 5, 8, 9 

0.214 

1,2, 3,4,5, 7,8, 9 

O 

1,2,3,4,5,7,8,9,10 

0 

1,2,3,4,5,7,8,9,10 


0.25 0.90 3 


1 

251.10 

2 

197.05 

3 

150.10 

4 

104.26 

5 

63.75 

6 

25.62 

7 

2.40 

8 

0.39 

9 

0 

lO 

0 


1 

1.4 

1 . 3.4 
1 , 3 , 4, 8 

1 . 3 . 4 . 8. 9 

1 . 3 . 4 . 5 . 8. 9 

1 . 2 . 3 . 4 . 5 . 8. 9 

1 . 2 . 3 . 4 . 5 . 7 . 8. 9 
1 , 2 , 3 , 4 , 5 , 7 , 8 , 9 , 10 
1 , 2 , 3 , 4 , 5 , 7 , 8 , 9,10 


0.25 0.95 3 1 


266.8 


4 


2 

3 

4 

5 

6 

7 

8 
9 

lO 


209.2 

159.34 

111.20 

68.13 

27.53 

3.05 

0.649 

0 

0 


1.4 

1.3.4 
1,3, 4, 8 

1 .3. 4.8.9 

1.3. 4. 5. 8. 9 

1.2. 3. 4.5.8. 9 

1 . 2 . 3 . 4 . 5 . 7 . 8. 9 

1.2.3.4.5.7.8.9.10 

1.2. 3. 4. 5. 7. 8. 9. 10 
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results corresponding to optimal design at network for groundwater 

quality monitoring. Here a refers to the degree of uncertainty 

for the corresponding CDFs of actual concentrations, ft refers to 

J 

the probability that the actual concentration will not exceed a 

particular predicted value at site j and time period k. P refers 

to the maximum number of wells permitted to be installed for 

groundwater quality monitoring. S. refers to the allowable limit 

of the pollutant concentration at site j. 

Figures 12 through 21 shows optimal observation well locations 

obtained as solution to the model for a = 0.1. ft = 0.75. S =3 

j i 

pci/1 far all j's with P varying from 1 to 10. Similarly Figures 

22 through 31 shows optimal well locations obtained as solution to 

the model for a = 0.25, ft = 0.75, S. = 3 V ci/l for all values of 

3 J 

j's with P vrying from 1 to 10. 

Increasing the reliability from 0.75 to 0.95 for a given 
degree of uncertainty (a), will result in a more conservative 
prediction of the spatial and temporal distribution of the 
pollutant concentration due to known sources. Figure 32 shows the 
optimal well location obtained as solution to the model for a = 
0.1, = 0.95, P = 2. Figure 33 similarly shows the optimal well 

location obtained as solution to the model for a = 0.25 ftl = 0.9, 

•r 

and P = 1 . 

These results clearly demonstrate that i) the degree of 

uncertainty <ot> in the CDF's of actual concentration and ii) the 

_k 

explicitly stated reliability of the predicted concentrations (3 
significantly affect the solution obtained from the optimization 
model. A measure of the degree of uncertainty is given by a 
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Sj =3 pci 
a = -1 
ft-- 0-75 

r-' 


■ Pollution source location 

© Potential observation well location 

(5) Observation well location according to 
optimal network design e 

Figure 12 Optimal observation well locations 
obtained as solution to the model 





Sj =3 pci II 
a = 0-1 
ft*- 0-75 

f =2 


■ Pollution source location 

® Potential observation well location 

(•) Observation well location according to 
optimal network design 

Figure 13 Optimal observation well locations 
obtained as solution to the model 
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Sj =3 pci /I 
DC =0-1 

£ k =0-75 
' J 


■ Pollution source location 

® Potential observation well location 

® Observation well location according to 
optimal network design 

V 

Figure 14 Optimal observation well locations 
obtained as solution to the model 


















Sj =3 pci II 


a =0*1 

A k = 0*75 
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f = 4 


■ Pollution source location 

® Potential observation well location 

® Observation well location according to 
optimal network design 

Figure 15 Optimal observation well locations 
obtained as solution to the model 
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Sj =3 pci / I 
a = 0-1 
/3 k = 0-75 

P =5 


■ Pollution source location 

® Potential observation well location 

® Observation well location according to 
optimal network design 

Figure 16 Optimal observation well locations 
obtained as solution to the model 
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Sj =3 pci / I 
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■ Pollution source location 

® Potential observation well location 

® Observation well location according to 
optimal network design 

Figure 17 Optimal observation well locations 
obtained as solution to the model 
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■ Pollution source location 

® Potential observation well location 

(•) Observation well location according to 
optimal network design 


Figure 18 Optimal observation well locations 
obtained as solution to the model 
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■ Pollution source location 

© Potential observation well location 

(•) Observation well location according to 
optimal network design 

Figure 19 Optimal observation well locations 
obtained as solution to the model 
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■ Pollution source location 

® Potential observation well location 

(•) Observation well location according to 
optimal network design 

Figure 20 Optimal observation well locations 
obtained as solution to the model 





Sj =3 pci 

oc = 0*1 

ft k =0-75 
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■ Pollution source location 

® Potential observation well location 

® Observation well location according to 
optimal network design 

Figure 21 Optimal observation well locations 
obtained as solution to the model 
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■ Pollution source location 

® Potential observation welt location 

® Observation well location according to 
optimal network design 

Figure 22 Optimal observation well locations 
obtained as solution to the model 

















Sj =3 pci 
oc =0-25 
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■ Pollution source location 

® Potential observation well location 

® Observation well location according to 
optimal network design 

Figure 23 Optimal observation well locations 
obtained as solution to the model 
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Figure 24 Optimal observation well locations 
obtained as solution to the model 
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Sj -3 pci /l 
a =0-25 

^= 0-75 
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■ Pollution source location 

© Potential observation well location 

© Observation well location according to 
optimal network design 

Figure 25 Optimal observation well locations 
obtained as solution to the model 
















■ Pollution source location 

® Potential observation well location 

® Observation well location according to 
optimal network design 

Figure 26 Optimal observation well locations 
obtained as solution to the model 
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Figure 27 Optimal observation well locations 
obtained as solution to the model 
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Figure 28 Optimal observation well locations 
obtained as solution to the model 
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Figure 29 Optimal observation well locations 
obtained as solution to the model 
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Figure 30 Optimal observation well locations 
obtained as solution to the model 
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Figure 31 Optimal observation well locations 
obtained as solution to the model 







































i 


I 



Sj =3 pci /I 
a = 0-1 
fiSo-95 

' J 

F =2 


■ Pollution source locotion 

© Potential observation well location 

(•) Observation well location according to 
optimal network design 

Figure 32 Optimal observation well locations 
obtained as solution to the model 
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Figure 33 Optimal observation well locatiops 
obtained as solution to the model 
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values. The explicitly specified reliability levels are given by 

|c- 

corresponding value of 

A higher value of a signifies a greater degree of uncertainty 
in estimating or simulating the spatial and temporal distribution 
of concent ratio ins. A higher value of signifies, a more 

conservative prediction of these concentrations, for a given 
degree of uncertainty given by a. 

A decision maker has to choose a particular value of <s, 
depending on the perceived uncertainty in the transport modeling. 
A specified value of a denotes an overall subjective 
quantification of the uncertainties in modeling the transport 
proc«.s. A choice of /I* on the otherhand, will reflect a «r« 
conservative prediction of concentration for an already perceived 
degree of uncertainty- values relate to the reliability of 

fc- 

prediction in the sense that a higher ft: value will ensure that 
the probability of exceeding a predicted concentration is lower. 

A higher value of a signifies a higher degree of uncertainty in 
the modeling of transport. It is however, possible to obtian 
more realistic estimates of a values to be specified by using 
conventional methods of uncertainty analysis such as a first order 
uncertainty analysis (Tung, 1986). 

The solution results demonstrate the feasibility of using 
this particular model linking a simulation and a chance 
constrained mixed integer programming aptimiztion model for 
optimal design of groundwater quality monitoring network. These 
results also elucidate the significance of using different 
uncertainty measures and reliability values. 
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CHAPTER 4 

SUMMARY AND CONCLUSION 


4 .. 1 SUMMARY 

A mathematical model linking a groundwater pollution 
transport simulation model and a chance constrained optimization 
model explicitly considering uncertaint ies of transport simulation 
is developed for designing an optimal groundwater quality 
monitoring network. To illustrate the performance of this model a 

T 

radioactive pollutant (Tritium, f-T 1 > is considered. Therefore, 
this developed methodology is particuvlarly applicable to 
radiactive waste disposal and groundwater quality monitoring. 
However, this method is easily applicable to other kinds of 
pollutants. 

The simulation model provides the information about the 
pollutant transport with respect to time and space. In this study 
the groundwater transport simulation model has been combined with 
a statistical generation technique in an effort to incorporate 
uncertainties in the model parameters. 

To provide a link between optimization and simulation model 
the response matrix approach has been utilized in which an 
external groundwater transport simulation model is used to develop 
unit response. An assemblage of the unit responses, a response 
matrix, is then used to simulate solute transport in groundwater. 
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for any given combination of pollution source location and 
magnitudes. Hence the problem of repeated simulation to provide 
input to optimization model is removed. In this approach, the 
response matrix itself is assumed to be uncertain. 

A chance constrained optimization model is developed for an 
optimal design of groundwater quality monitoring network. The 
role of the chance constraints in the model is to introduce same 
degree of reliability in the predicted values of actual 
concentrations at specified locations. The optimization model 
incorporates the cumulative distribution function <CDF> of the 
actual concentration far each potential groundwater quality 
monitoring well location and for each specified management time 
period. The CDF's can be obtained by randomly varying the 
response matrix, that represents the response of the subsurface 
saturated zone for a given input of parameters. 

The main objective of the chance constrained optimization 
model is to maximize the probability of detecting contamination at 
locations, where the standard is exceeded with due weightage given 
to its degree of exceedences. The exceedence values are computed 
in terms of the predicted values of concentration expressed as a 
function of specified reliabilities. The optimal design of 
monitoring well network restricts the total number of monitoring 
wells to a presepecif ied upper limit and at the same time takes 
into account the uncertainties of predicting salute transport in 
groundwater . 



The performance of the developed model was evaluated for a 
number of scenarios with different degrees of uncertai inties and 
different specified reliability values. These evaluations 
demonstrate the appl icability of this model for design of 
groundwater quality monitoring networks, with inherent 
uncertaint ies in transport modeling. 

4.2 ENGINEERING APPLICATION OF THE DEVELOPED METHODOLOGY 

The methodology developed in this study is primarily useful 
for monitoring the temporal and spatial movement of a contaminant 
plume in groundwater. This detection is the first step in proper 
management of groundwater contamination. The optimal monitoring 
network ensures that the monitoring wells are installed at 
locations where these wells are able to detect contamination with 
high probabilities. Uncertaint ies present in the modeling of the 
contaminant plume movement require that the monitoring design 
network incorporate probabi lities or reliabilities of detection. 
Given the upper limit on the number of monitoring wells to be 
installed generally based on economic considerat ions, the optimal 
design locates those sites where the resulting concentrations are 
probable to be very high. Thus it ensures an efficient network 
for monitoring while considering the uncertainties in parameter 
estimation, measurement errors, and other modeling uncertaint ies* 

Depending on the actual purposes of groundwater withdrawal, 
such as domestic water supply, irrigation water supply, industrial 
uses and others. This model is capable of using different 

allowable standards for contaminant concentrations at different 


locat ions . 


Therefore, if for example, groundwater at a location 
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is being used far drinking 

^ P^rp OS8 s . i 

. . . . . ’ ^wer value of allowable or 

standard concentration Car> ' 

wore efficient monitoring location, to ensure 

. . . ^ *° c ation. Also, different 

weights can be introduced an 

... he object ive function, to reflect 

the importance given to the 

^sed.ncs of a giv6n standard at a 

given location. 

In particular, this . 

roethod Ql win . . . . 

yy il be useful m monitoring 

the effects of radioactive Da n 

utant disposals from nuclear power 

plants. Also, the opti mal 

i to ring network can be used to 

detect and subsequently raan 

contamination of the groundwater 

due to accidental disp Qsal , 

radioactive wastes. 

4. 3 CONCLUSIONS 

1. A mathematical model xint ■ 

109 a groundwater solute transport 

simulation model and a rh 

a ^ce constrained optimization model 

is developed. The optimi, . . 

-a -ion model explicitly considers 

uncertainties of t r , 

r simulation and incorporat 

specified reliability 

of Predicting the 
concentration in time and space. 

The model is evaluated f QP .. 

ifferent degree of uncer taintie 

in the modeling process l. 

by introducing CDF's of actua 
concentrations as constraint- in 

ta in the optimization model. 

The solution results show th*+ *-h 

the optimal well locations ar< 

significantly dependent Qn w, ^ 

Qn the decree of 

the -specified reliabiliK, 

^ values. 


es 
actual 


*L * 


3 , 


uncertainties and 
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APPENDIX - lA] 


FORTRAN Programme for 
Statistical perturbation of the 
[R] matrix. 



C £ T /a T 1 *:■ f i € A l PEfcTUkBAIION OF THE CRT HATRIX 

i i i 

C a| * 4 4 * t| * ! * k * 4' * * % 4 A ) 4 * » 4 * 4 .* * » j|. 4 . * * 4 4 4 4 *1 * * ♦ * ♦ 4 *f= * 

DIME NS ION i t 1 0 ) , C < 1 0 UO ) 

DIMENSION M 1000, Ho . CN < t 0 0 0 , 1 0 ) , CN 1 ( 1 000, 10) 
OPEN t UN IT -SO , F IL E * ' dh i men . pm ' > 

OF £ N < UN I T » 3 1 , F IL L« ' sun da t ' 1 
OPEN (UN I Ucl,FU£='hc") 

OPE N l UNIT = ££ .. F ! L £ - ' K3 ' ) 

OP t N < UN-1 T * 3 0 , F I L E •- 'HI ' ) 

OP t N ( u N IT - c 3 . F I l E ~ ' N 4 ' ) 

OPEN 1 UNIT * , FILE* * HE ' ) 

OPEN 1 UNI 7 *3E , F i L£* ' H£ ‘ ) 

OPEN t UNI T*3£ , F I L £ = ' H 7 ’ ) 

OPEN t UN 17*37, FILE*' H8‘ ) 

OPEN ( Un I T = 38 ,FILE* ' HS* ' ) 

OPE N ( UN 1 1 * 3 V , F I L .£'* ; H I 0 ‘ i 

ALPHA* 0 £E 

RHEA N = 0 , C 

IP* ! 7 9 S 4 

DO 1*1.10 0 0 

RE AD ( SO , * ) ( h ( 1,3) ,3 = 1 , 8 ) 

ENDDO 

R E A D < 3 1 , *i > ( 5 ' I ) , 1*1 , & ) 

DO e 0 0 M * 1 , 1 C* V 
DO 1*1,1000 
DO J - 1 , £ 

& I & M A *■ R <1,3 1* ALPHA 
CALL SPANDt'IP) 

15 R1*RAND() 

R 2 * R A N b ( ) 

Y * l - c . * At OG i R 1 ) ) * *0 . E* COE < $ . £83*R£ ) 

CNl I ,3 ) "SIGMA* Y+RMEAN 
cni a , j )*R( i , j s-rcm 1 , J> 

IFlCNHI.J) .LT.0.0) 00 TO IE 

IP * IP+i 00 

ENDDO 

ENDDO 

DO I *1,1000 

cm = 0.0 

DO £0 J*1 ,& 

C ( I ) * C ( I ) +CN 1 ( I , J ) * S ( J ) 

20 CONTINUE 

ENDDO 

URITE(30,*)<C(I >, 1=10,100,9) 

WRITE i 21 , * M C( I 5 , 1*1 1 0,£00, 9) 

UR I TE <££,*) (0( 1 ) , 1*21 0,300,9 ) 

WRIT E( £3, *)<C( I ), 1=310,400,9) 

UR I TEC 34, ■» ) (Cm , 1*41 0,500, 9) 
WRITE(3&,*)(C(I), 1=510,600,9) 

WRITE (3b , * ) i Cm , 1*61 0,70 0 , 9 ) 

WRITE ( 37, * ) (C( I > , 1=71 0,800,9) 

WR I TE( 38, *) (C<I) , 1-810,900,9) 

WRITE ( 39 , * ) (Ct I ) , 1=91 0, 1 000,9) 

IP-IP+ 1 00 
200 CONTINUE 
STOP 
END 



